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Abstract. The classic theory of electron field emission from a cold metal surface 
due to Fowler and Nordheim (FN) is re-examined and found to violate the validity 
criteria for the WKB approximation, for electric fields greater than about 1 V/fim. 
In this study we shall examine the complete solution without invoking the WKB 
approximation in order to assess the reliability of the FN theory as widely used for 
the interpretation of experimental data. Particular problems occur when the barrier 
height (and therefore indirectly also the width) is significantly reduced by the image or 
some external effects. Further refinement of the theory will be discussed by considering 
the effects of screening, which can be one mechanism for the barrier height reduction, 
in addition to the widely known negative affinity of diamond like carbon systems. A 
comparison with experimental data from carbon field emitters shows that the enhanced 
current found in this paper may provide an explanation for strong field emission 
observed recently in carbon based samples. 
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Introduction 

The classic theory of electron field emission from a cold metal surface due to Fowler 
and Nordheim (FN) QjJ |2l El was one of the earliest applications of wave mechanics 
since its foundation in the early 1920's. With increased interest in field emission 
processes, especially for use in display screens it is worth examining closely the 
foundations of this theory. As well as metal-vacuum interfaces, proposed applications 
of the theory includes graphite-vacuum, graphite-insulator-vacuum, carbon-nanotube- 
vacuum and Schottky junctions Since our approach can be generalized to other 
cases of interfaces between dissimilar media, it may have application to the systems 
discussed by Mott and Gurney |0] and Herring |7|. These other cases might include 
vacuum breakdown and electrical forming processes jH] . In section [S] we show that our 
new analysis removes discrepancies of many orders of magnitude between theory and 
experiment. 

In the FN theory, the role of the Fermi surface, Fermi-Dirac statistics and in 
particular the unique property of quantum mechanical tunnelling were key concepts. 
The essentially equilibrium distribution of electrons is modeled by a one dimensional (x 
dependent only) Hamiltonian: 

H = ^ + V(x), (1) 

where all energies are measured with the same reference as the effective potential V(x). 
Then the current density J of the emitted electrons is given by an integral over the 
energy W: 

/oo 
P(W)dW, (2) 
-W a 

where P(W)dW gives the number of electrons within dW that emerge form the 
metal per second per unit area, with W a an appropriate lower energy reference. The 
quantity P(W) is easily expressed as the product of a quantum mechanical tunnelling 
transmission coefficient T(W) and a Fermi type supply function N(W): 



N(W) = — ^— In 1 1 + exp |-(— ^)j j , (3) 



W-e F . 



T(W) = exp 



(4) 



in which ep is the Fermi energy. In view of the complexity of the Schrodinger equation, 
the transmission problem was originally treated within a WKB approximation or 
variants of it 9\: 

{^(V(x)-W)j d 

in which X\ and x<i are the classical turning points of the potential V(x). It is a well 
known result however that the validity of the WKB approximation involves a criterion, 
see section^ that can be derived from the requirement that the variation of the deBroglie 
wavelength of the electron must be less than the dimension of the region of tunnelling 
[HI] . This restriction on the validity of the FN theory was unfortunately not discussed in 
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many of the earlier papers, including the classic paper of [3]. The limitations of the WKB 
approach are well known, however, in the treatment of thermionic emission; a common 
technique for improving the calculation [IT] is to make a parabolic approximation to the 
top of of the barrier. Nevertheless, FN theory has been applied in regimes where it is not 
valid. In section El we shall examine the image- free problem exactly whose solution in 
one dimension is amenable to analytical treatment in terms of Airy type functions. Some 
preliminary assessment of the FN- WKB can already be made for this case, which is a 
preamble to the full image problem which we shall treat in section E]by various analytical 
methods. Interestingly the tunnelling problem close to the barrier top appears to be 
amenable to analytical treatment as well and has not been published to the best of our 
knowledge. We shall note that the failure of the WKB criteria coincides quite closely to 
this regime, hence the availability of analytical results are extremely welcomed. Next 
we shall note that the role of screening has also been omitted in the original FN theory. 
We shall introduce the additional screening potentials in section |U and thereby propose 
that this is one mechanism for a significant barrier reduction, in addition to negative 
electron affinity attributed to diamond like carbon systems. In section El we shall apply 
our theory to the analysis of experimental data for carbon field emitters |12j. Regardless 
of reasonable assumptions of point tip field enhancement effects, we shall see that the 
FN theory does not even get close to the order of magnitude agreement with the data 
by comparison with ours. We shall conclude in section |H1 with a summary of our results, 
postponing further solutions of the problem involving screening to a future paper. 



1. WKB Validity criteria 

Following it is straightforward to derive the criterion for the validity of WKB: 

p 3 » mheF, (5) 

where F is the external field and we shall ignore the electron image for the moment. 
Here p is the electron momentum so that: 

(2m((f) - eFx)) 3/2 » mheF, (6) 

where <fi is the work function. Assuming a carrier effective mass equal to the free electron 
mass and inserting numerical values appropriate to a field in V/fim and a work function 
in eV, (which shall be our units throughout this paper) gives 

C(F, <P) = 0-046^ «1, (7) 

in which the tunnelling length x has been set to zero, representing the best case for 
WKB, the very top of the barrier (x is typically of the order of nm). Typical numerical 
values for ( are given table ^ where we see that the WKB approximation begins to 

| Since image effects are less important for silicon based Schottky junctions due to the large dielectric 
constant e s » 11.9, these results are actually very relevant to the latter, which however we shall not 
take up here in this paper. 
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Table 1. Values of the function <f). For WKB theory to be valid, ((F, (j)) « 1. 



F(V/nm) <j) = 1 eV <j> = 5 eV = 10 eV 



1 


0.046 


0.021 


0.015 


10 


0.099 


0.044 


0.031 


100 


0.213 


0.096 


0.068 


1000 


0.460 


0.206 


0.146 



be questionable for fields above about 10V//im for typical work function values and is 
certainly seriously questionable in the range close to dielectric breakdown fields which 
are typically 10 3 V/fim . A smaller effective mass such as for p-Si where (m* = 0.6 for 
holes) would increase £ further. Later on in the paper we shall show further anomalies 
with the FN theory when we are close to this limit of validity. In order to proceed and 
be able to assess the approximation further we shall have to resort to a full solution of 
the Schrodinger equation which we shall begin in the next section by first ignoring the 
image potential. 

2. Exact solution without electron images 

We shall start with the image free tunnelling problem given by the Schrodinger equation: 

tP"(x) + -z[W -V(x)}ijj(x) = 0, (8) 
n 

where W is the energy eigenvalue of the Hamiltonian Hq with the potential given by 
0: 

f-W. if, <0 -region 1; 

I V — eFx if x > - region 2. 

Here V) under normal circumstances is zero which yields the work function as the 
barrier height and F is as before the applied external field. For later convenience we 
shall recast Eq(JHJ) into the more compact form: 

ip"(x) + (e + ax)i/>{x) = 0, (10) 



with 



a 



2meF 

IT 

2m(W - V 



F x 2.55989 x 10 7 (/W) 



-3 



. - 2 ~~ (W - V ) x 2.55989 x lOVm)" 2 , (11) 

in which F is measured in V/ /zm. To the best of the authors' knowledge this elementary 
quantum exercise does not seem to have been fully documented in the literature, 
although variants of it have been scattered around O EU H3 QSj. Following Landau 
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and Lifshitz ^U|, we shall transform to the new variable £ = (x + e/a)a 1//3 . Then 
equation (fTU|) in region 2 becomes the typical Airy's equation §: 

<(0 + ^(0 = o. (12) 

However some care needs to be exercised in selecting the appropriate type of Airy 
functions, since the standard Ai function has asymptotic form [TU1 117j: 

m = Ai(-0 « T 1/4 sin(^ 3/2 + tt/4). (13) 

As this also contains a reflected wave, it is not appropriate to our boundary conditions. 
The correct approach is to revert to the original 1/3 fractional Bessel functions definition 
for the Airy functions and then construct appropriate linear combinations to yield the 
forward travelling wave solution only. This approach of course leads to a modified 
appropriate Hankel function which as we shall see has the correct asymptotic form 
required. It is this construction that seems to be rather rarely seen in the literature 
except perhaps in the WKB derivation of Moll [TB] . With this proviso, we can write 
down the solution for the wave functions of equation (fT2j) as: 



e lkx + re tkx in region 1; 

' ' - e- i7r/3 w 2 (0] in region 2, ! 1 ' ' 

where k = (2m(W + W a ) /h 2 ) 1 ^ 2 and the u\ and 112 are given by Bessel functions of order 
1/3, here defined as: 

«i(0 = ^K) v V- 1/3 (^ 3/2 ), (is) 

and 

= \{^) 1/2 J^{\e'\ (i6) 

respectively. Note however that analytical continuation to negative £ requires some care: 
see Appendix A. It is then straightforward to obtain the reflection and transmission 
amplitudes by matching the wave functions and their derivatives at the interface x = 0. 
We provide the results here in terms of i?i/ 3 = u% — e^ m ^u 2 which is a slightly modified 
Hankel function 

2 

t = 



tfl/3(£o)-^0 1/3 #W£o) 



1/3 * 

H 1/ ^ )+ i k 1 /3 H[ / .^ ) 

H 1/3 ^ )-iki /3 H' 1/3 ^ y 

where £0 = £(0) = e = e/a 2//3 while k^ 3 = a 1//3 jk. An important check must be made 
to these formulas which requires that the transmission coefficient: 
1 1/3 

T = ^ r |t| 2 , (18) 

§ The advantage of using the Airy function formulation is that it naturally takes care of the long range 
boundary condition posed by the electric field. 
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J- asymp ~ t2/3i- ^ ' 1 ' 



and reflection coefficient: -R = |r| 2 satisfy the unitarity property T + R = 1. This is 
easily checked to be the case when appropriate use is made of the Wronskian identity: 

um' 2 - u 2 u 1 = — (19) 

For sufficiently small fields, we can obtain an asymptotic expansion (where £o — > oo and 
fco — ► 0) that yields: 

:i + ^o 2/3 i£oi) 

the prefactor being a new result not obtainable within WKB. For convenience we set 
W a = 20 which for of order 1 eV is adequate as a lower bound (at room temperatures), 
the exception being very narrow band semiconductors which we are not concerned with 
here. We shall at the moment ignore any possible barrier height lowering (i.e. image 
free case) so that Vq = is the vacuum level and consider tunnelling near the Fermi 
level, so that W = — and thus fc(^ 3 |£o| = 1 . Then we see that the above prefactor 
differs from the WKB result by a mere factor of two i.e.: 

Tasymp ~ 2 6 3 . (21) 

That this asymptotic result is valid for only small fields is clearly seen now for otherwise 
if we naively allow the field to be large, then we could arrive at an unphysical result 
in which T > 1 which is absurd. In fact this occurs for 0=1 eV when F reaches the 
order of 9.9 kV/^im. Hence the use of the WKB formula for large transmissions without 
knowledge of the prefactor in the transmission can be dangerous. Unfortunately this 
caution does not seem to have been exercised by Good and Miiller nor by Murphy 
and Good JB]- The latter have derived widely used criteria JH] for the validity of FN 
theory on the basis of the WKB integrand properties equation (@J) alone, which is thus 
a logically inadequate procedure. If we examine the classic FN formula for the emission 
current: 

>- , t ^:T m ^-(g^W >< a ~- m 

in which t(x) and v(x) are given in terms of elliptic integrals, for = 2 eV say, the 
exponent practically vanishes at F m 3000V//zm. Hence the entire emission current 
comes from the supply function which is the F 2 factor and, as discussed above, this can 
be untenable. In fact it is well known from elementary quantum mechanics that the 
transmission coefficient is in general not unity at the top of the barrier [2*U] . 

Later on, in sectional we shall see that an exact solution is obtainable for the full 
image theory at the top of the barrier, in which the transmission coefficient has a finite 
value less than unity, but nevertheless with significant tunnelling current contributions. 
As far as we are aware, this solution has not been noted in the field emission literature for 
the last 40 years. In figure ^ we provide two curves that compare the exact transmission 
coefficient equation (|18j) for the image free case against the asymptotic formula equation 
(12 lj) close to the Fermi level for very large fields where the latter becomes problematic. 
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Figure 1. Log-linear plot of the exact transmission coefficient (solid line) equation 
(|18|l vs the asymptotic formula equation (|21|l (dashed line) close to the Fermi-level in 
the image- free case for large fields. F is in units of V//im. 



Fortunately these fields are very large and being much greater than typical breakdown 
fields this failure is unimportant in practice. We note that for smaller fields figure El the 
WKB formula is quite good, in spite of its failure to satisfy our earlier criterion of section 
^ The accuracy of the WKB for small fields is however due to us staying away from 
the top of the barrier. We now compare equation ()21|) with the exact result equation 
(jl8|) as we approach the barrier, i.e. Vq —* —<p. One mechanism for this to occur could 
be due to a large image term i.e. Schottky effect, but there are other screening effects 
as well as we shall see later. For simplicity we shall not include the image effect in the 
Schrodinger equation yet but merely consider the effect of approaching the barrier top 
in the image free case first. Here we can compare the results of three formulas, the exact 
equation (|TH|) . the asymptotic formula equation (|2"UJ) and the standard WKB result in 
which the prefactor in equation (J2"U|) is missing: these are shown in figure El where we 
have set Vq = —0.95 <p. Here we see that for quite small fields F < 1001^/ /im, significant 
departures from WKB are found. These results motivate us to further include the 
Schottky image term in the Schrodinger equation which we take up in the next section. 



3. Exact solution including electron images 



We now turn to the quantum tunnelling problem including the effect of the image term 
as in the original FN theory. Then equation (fTU|) becomes: 



ijj"{x) + 



e + ax + 



= o, 



(23) 



where the image potential is characterized by the parameter (3 = me 2 /(2h 2 
Transforming to the parameter £ as before we now have the equation: 

a 



r(0 + £ + 



^(0 = o, 



(24) 
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Figure 3. Log-linear plot of the exact transmission coefficient (solid line) equation 
(|f 8|l vs the asymptotic formula (dashed line) equation (|20|l and vs the standard WKB 
formula (dotted line) close to the top of the barrier (Vo = —0.95 <f>) in the image-free 
case. F is in units of V//zm. 



where a = /3a 1 ^ 3 , and e = ea 2 / 3 . Clearly the image-free case a = reduces to the 
usual Airy differential equation equation (|T2*j) and the above can be rewritten as: 

(C - e)V"(0 + [m - e) + = 0. (25) 

We shall in fact further rewrite this as: 

r(e-Ai)(£-A 2 )- 



V>(£) = o. 



(£ - e) 

in which the two roots of the quadratic are given by: 



Ai 



1 



-[e + Vei-Aa] and A 2 = -[e - V? - 4a]. 



(26) 



(27) 



Let us first discuss the WKB approximation ^U], which is obtained by approximating 
the wavefunction as: 



p(£) 



cxp 



Z7T 

T 



(28) 
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where: 

rtfl=[ «-*y> r (29) 

and the lower limit in the integral in equation ()28|) is the smaller of Ai or A2 which 
depends on the sign of e i.e. above or below the barrier respectively. The evaluation 
of the integral in equation (J2"%|) is then a straightforward exercise in Elliptic integral 
reduction which we need not repeat here El EE]- Equation (ffijj) itself does not 
appear to have a closed form solution in terms of known functions although several 
methods of solution are available depending on the parameters. We postpone these to 
the appendices as we do not intend to pursue detailed calculations here, but merely to 
demonstrate the regimes where WKB is in need of correction. For this purpose we can 
concentrate on the region close to the top of the barrier at £ = a 1 / 2 + e. Let us start 
at the barrier top itself for which WKB predicts T = 1 since in this case v(x) = in 
equation (|22|) . Then equation ()26|) simplifies to the form: 



<(o+[ ( ^5| 2 ]^) = °- ( 3 °) 

In view of the magnitude of the parameter e it suffices to obtain a solution in the limit 
when this parameter is large. For = 1 eV and 100 < F < 1000 V//xm, e ranges from 
3 to 14. In particular for the matching conditions we are only interested in the solution 
near to £0 = £ which corresponds to the origin x = so that we can transform to the 
variable: y = £ — e which is small for our purpose y << e so that equation (JHOJ) can be 
approximated by: 

#"(2/) + (y + e/2) 2 ^(y) « yip"{y) + mf^iv) = o- (31) 

The appropriate solution for the above is then given by the Hankel function of order 1 
which upon transforming back to the original variables £ is given by: 

HO = e W4 /| - ^H^m ~ ^f' 2 )- (32) 

We can now obtain the transmission coefficient as was done in the previous section by 
matching the wave function and its derivative at the origin. The result is now given by: 

2e 2 kl /3 

T ~ l^o)| 2 + e 2 ^fc 2/ V(£o)| 2 (33) 
where the derivative of ip is easily shown to be given in terms of the Hankel function of 
order 0: 

^'(£) = WJ e^^e-e) 1 /*) (34) 



As usual we can verify that this solution satisfies unitarity T + R — 1 through the use 
of the appropriate Wronskian identities for the Hankel functions. That this is satisfied 
exactly in spite of our approximation in equation (}3T|) means that we have picked out 
the essential features in the theory. Equation gives zero transmission at the top of 
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the barrier because of a weak logarithmic singularity (typical of the 1/x potential near 
the origin) in ip'. We may write 

„ 4:71 



~ e 3 / 2 (ln(5 1 /2e))2 ( 35 ) 

where 5 is a renormalization constant, which we can arbitrarily set to 5 << e. Physically, 
we may think of this as correcting for the situation in which the electron and the image 
charge are both at the interface by limiting the distance of closest approach to atomic 
dimensions. Equation (|35|) corrects for the error made by the WKB approximation at 
the top of the barrier. As a convenient choice of 5, we require that it is chosen to limit 
the closest approach to the origin as x = a = 1A || so that by the definition of £: 

10 _4 F 

S = ——. (36) 



Before proceeding further we shall briefly mention the results when we move slightly 
away from the top of the barrier. In this case the two roots become: 

Ai « 1(1 + 0-) and A 2 w|(l-«r) (37) 

where a = (1 — 4a/e 2 ) 1 / 2 and in the limit e large reduces equation (|3~Tj) to: 

y^"{y)+^{y) = 0. (38) 

where 7 = ±e(l - a 2 ) 1 / 2 . This poses no difficulties and the method outlined above can 
be adapted to the solution. We summarize the appropriate formulas without further 
ado: 

8l 2 kl /3 



(39) 



where: 



and 



m = e^ /4 /| 2 7 (£ - e) 1 / 2 J ffW ( 2 7 (e - e) 1 / 2 ) (40) 
i>\0 = e m/ \[l (2 7 ) 2 ^ 1) (2 7 (e-e) 1/2 ) (41) 



In figure 0] we compare the transmission coefficient equation (J33|) versus the FN 
formula equation (|22jl less the prefactors. We see that the FN formula is not only in 
error as we trespass the barrier giving a non physical value of T > 1, for the region 
near the top of the barrier (around F = 700\^//zm here), it drops exponentially fast as 
opposed to the exact solution. 

Using equation (|39j) we can estimate the current near the barrier top which occurs 
for F ~ 700 V//x in this case, remembering that we have considered 0=1 eV and 

|| By way of comparison, the typical tunnelling length C — f aus i n the range of 10 4 A to 10A for 
1 < F < 10 3 V/fjxa. 
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Figure 4. Log-linear plot of the exact transmission coefficient equation (|33|l f solid 
line) vs the FN formula (dashed line) i.e. equation less the prefactors in the full 
image case. F is in units of V/fim. 



W a = 2(f). Hence for temperatures much less than ep/ks, the current is obtained from 
the integral: 



which has to be evaluated numerically. We can however make an estimate using the 
approximate expression equation (|35|) . We can set the Fermi level at the top of the 
barrier and integrate only over a small region of order S near to it, since outside this 
region the contributions are exponentially small. We choose 5 according to equation 
(jnnj) as before and thus the e _3//2 pre-factor dominates in the transmission coefficient 
since the logarithmic term is slowly varying in the region of integration. As we take 
W a /cf) = 2 which is large compared with 5, we now arrive at the expression for the 
current density as: 



Note that the field dependence of the current is in this case F 3 which differs from the FN 
prediction of F 2 and could be experimentally discernable in view of its non-exponential 
dependence. We can now compare equation ()43|) with the standard FN formula equation 
(J221) see figure El 

Having shown that we can have a large (non-exponential) current in spite of a 
smaller field than in the image- free case, due to proximity to the barrier top, we shall next 
examine other mechanisms by which the barrier lowering can occur, whose magnitude 
is comparable with, or might exceed the Schottky image effect. The effect to consider 
is due to the dielectric screening of the applied field F which we shall take up next. 

4. The Effects of Screening 




(42) 



J « 0.1674 <p~ 3/2 F 3 A cm" 2 



(43) 



Screening is an important property of metals or dielectrics, and the mechanism itself 
plays a crucial role in the previous section in the form of the image potential which is 
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Figure 5. Log-linear plot of the tunnelling current near the top of the barrier (which 
occurs for fields of order 700 V//im here) (solid line) equation (|43J) compared with the 
standard FN formula (dashed line) close to the top of the barrier in the full image case 
equation J22J). F is in units of V///m. 



a screening effect. In recent years there has been a revival of interest in screening, [21] 
especially in semiconductors, where surface effects can significantly alter the properties 
of devices. On the surface the screening length, typically of the order of Xtf from 
Thomas-Fermi theory ranges between the order of a few A for metals to a few /zm for 
semiconductors or even higher to tenths of mm due to surface effects since it is the charge 
carrier density at the surface layer that matters. The (Laplace) screening problem for a 
flat surface is easily solved for an external field F . This modifies the external potential 
V(x) in equation Q by having a non-zero V such that V = —Xtf^F. We shall neglect 
an exponential penetration field inside which in fact modifies the inner potential to 

V{x)= -W a - X TF eFe x/XTF for x < - region 1 

= — Xtf^F — eFx for x > — region 2 (44) 

The latter leads to a 2D electron gas which requires a self-consistent theory to treat 
since this in turn affects our Xtf as already mentioned. Further the proper inclusion 
of screening must now include its effects on the image potentials which pick up Yukawa 
type terms. It suffices here to estimate the barrier lowering from our expression for Vq. 
Since F is in V//zm then a Xtf ~ 1 f^ m w iU completely wipe out a barrier of 1 eV. For 
3D this corresponds to a carrier density of order 10 25 m -3 but for a 2D layer this can be 
down by a 2/3 power i.e. 10 17 m -2 . If we include the effects due to a diminished electron 
affinity as in the case of diamond-like carbon bonds, then we see that the present theory 
provides a strong candidate as an explanation for the high current densities in carbon 
granular samples recently observed in materials for field emission displays 4J, as we 
shall see in the next section. 
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Figure 6. Log-linear plot of the tunnelling current vs electric field F, assuming that 
a significant barrier lowering mechanism (Vo = 0.7(f)), such as dielectric screening (see 
text) has occurred. The solid line is our theory, the dots are experimental points from 
Tuck et al ^2] while the FN theory is given by the dash line. In addition a large field 
enhancement factor 7 = 10 has to be assumed for the FN formula equation Our 
result compares much more favourably with the experimental points than does the FN 
theory. 



5. Relation to Experiments 

Comparison with experiments such as for carbon field emitters ^2] can be made by 
assuming that there is a barrier lowering mechanism such as the dielectric screening as 
mentioned above. The tunnelling current will then be assumed to be given by equation 
(J43|) but we shall include a factor a to account for granularity or inefficient screening 
(and hence barrier lowering) which diminishes the current. This modifies equation (}43|) 
trivially to: 

J « 0.1674 0" 3/2 a 3 F 3 A cnT 2 (45) 

In figure El we have plotted the results which shows that for a granularity factor of 
a = 0.15 our results compare favorably with experiments ^2] bearing in mind that our 
results are only estimates and not accurate numerical evaluations^. In particular the 
parameters as given in the figures have not been fine tuned in anyway. On the contrary 
without assuming a significant field enhancement due to pointed tips, the FN theory is 
out by numerous orders of magnitude. However even with the assumption of a large field 
enhancement factor 7 = 10 which is reasonable for ellipsoids [23 J, and a barrier lowering 
Vq = —0.7(f), the FN curve figure El does not in any way resemble the experimental data, 
a feature we have found to occur quite commonly in analyzing data from the carbon 
field emitting samples. Note the vast order of magnitude differences and the fast drop 
in the slope of the FN curve which is shown more clearly in the inset to figure El These 
results clearly show that our theory has captured some essential physics associated to 
tunnelling near the barrier top which may be relevant to the carbon field emitters. 

% Finite temperature effects should also be included in a more accurate calculation, providing essential 
improvements to thermionic field emission theory |T51 1 19| . 
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6. Conclusion 

In conclusion we have presented a closer look at the traditional FN theory of field 
emission and found that the latter is untenable for calculations of emission currents 
when significant barrier lowering occurs. A full quantum mechanical solution of the 
problem provides us with an estimate which, coupled with a proposed barrier lowering 
due to screening or negative electron affinity, constitutes a reasonable order of magnitude 
estimate for the high current densities in samples with embedded particles, not available 
via FN theory. Further work is necessary to provide more quantitative comparisons 
noting that the granularity of the systems would lead to less total current density 
although on the other hand some field enhancements due to shape effects 123] not 
considered here might partially compensate for this loss. 
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Appendix A. Analytical continuation of Bessel Functions 

Confusion arises with Bessel functions of order 1/3 due to special care needed to treat 
the phase factors: some popular mathematical software uses incorrect continuation 
formulae. Firstly our convention is the same as Watson's so that the appropriate 
continuation formula is given by: 

_ j e- iU7r/2 J u (ze i7T/2 ) if -7T < arg(z) < vr/2 

lv{Z) = { e^l 2 J v {ze-^' 2 ) if vr/2 < arg(^) < vr. (A-1) 

An alternative approach is to perform the continuation for negative z as: 



3 , 



y/n sjx J_i/ 3 (|x2) i r— -2i 

«1 = 3 " = 3 V7T V M J -V3(^- 

o I I - 



= -e"/ 6 A^|L 1/3 (^), (A.2) 

which is incorrect. The negative argument prevents the use of the first form of 
equation (|A.l|) . Instead the branch cut properties requires the continuation formula 

Eh 

J u (e im7r/2 z) = e im7r/2 J u (z) (A.3) 
which holds similarly for the I u functions. Then: 

u i = g = 3 Vtt y\x\ J~i/3{— \x\*) = -^Tr^\x\I„ 1/3 {— — ), (A.4) 

without any complex phase factors. Similarly U2 picks up a minus sign: 

u 2 = 3 = 3 Vtt yj\x\ J i/3{— Fl a ) = 3 - Vtt v|x|/i/ 3 (— 3— )■ (A.5) 

Hence one has to manually continue the J v to the I v functions for negative arguments 
when if using the incorrect continuation, equation (jA.2|) . We have found that 
Mathematica v4.1 uses equation ()A.2j) . 

Appendix B. Power series solution for the full image problem I 

The generic form of the full image problem is given by the following second order 
differential equation 

if" + (e + ax + ^)ip = 0. (B.l) 

As before a straightforward linear transformation: £ = [x + —)oc* leads to the following 
form: 

QJ 

r + ^ + ^—^ = 0, (B.2) 
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where a = f3a~^ , and e = ect -2 / 3 . Clearly the image free case a = reduces to the 
usual Airy differential equation. Equation (jB.2|) can be rewritten as: 

(£-e)< + £(£-e)^ + c^ = 0. (B.3) 

Unfortunately this equation does not belong to the standard hypergeometric class and 
cannot therefore be treated by the usual mathematical physics functions. However using 
the standard Frobenius method one can obtain the following power series solutions. Let 



u 



(0 = E *nC +n , (B.4) 



m 

n=0 



then equation ()B.4j) leads to the following 5 term recursion relation for the determination 
of the coefficients a n : 

a n _i(m + n — l)(m + n — 2) — a n e(m + n)(m + n — 1) 

+ a n _ 4 - ea„_ 3 + «a n _ 2 = 0. (B.5) 

All the coefficients are well determined in terms of the coefficient a . The first term 
satisfies 

em(m — l)a = 0, (B.6) 

giving m = or m = 1 which shows as expected that we shall have two linearly 
independent solutions. Carrying on we have a x = and 

a 2 = —-, r • (B.7) 

e(m + 2)(m + l) V ; 

Notice that a± = regardless but a<i is only finite for a^O, i.e. (3 ^ 0. The following 
coefficients can be sequentially determined: 

_ a 2 (m + 1) a 

3 e(m + 3) (m + 3)(m + 2)' 1 ; 

Q _ a 3 (m + 2) | oq | a «2 g ^ 

° 4 ~ e(m + 4) e(m + 4)(m + 3) e(m + 4)(m + 3)' 
and so on. The general 5 term recursion relation can be re-written as 

1 



(m + n){m + n — 1) 

1 all 
-a n _i(m + n — l)(m + n — 2) + — a n _ 2 - «n-3 + -a-n-4 ■ (B.10) 
e e e J 

Unfortunately there is no way to simplify this any further. When a = (image free 

case) the above in fact reduces to a two term relation given by: 

a n = —- n 3 r. (B-ll) 

(m + n)(m + n — 1) 

This immediately leads us to the two independent solutions, for m = 
= x + g (-D-1.4 7 (3,-2) 

= 5 ^ 1/a ^(|M. (B.12) 
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and for m = 1 

(e) = < f (-ir2.5.8...(3,-l) 
llU ^„ti (3n + l)! 4 

= 3 1 / 3 r(4/3)e 1/2 J 1/3 (^ 3/2 ), (B.13) 

which are the series expansions for the Bessel functions of order 1/3. The oscillatory 
nature of these functions means that the convergence is in general poor and some more 
sophisticated re-summation method must be used to obtain in particular the asymptotic 
properties of these functions. In the next appendix we shall show just one method as 
to how this can be achieved. 

Appendix C. Power series solution for the full image problem II 

An alternative solution for the full image problem which is an effective re-summation of 
the previous power series solution that will allow us to perform numerical calculations, 
will now be shown. The approach is to first use a Laplace transform on equation (jB.ljl . 
Let 

m = I <Pf(t)dt, (C.l) 
Jc 

over some contour C, then the resulting differential equation in ip becomes converted 
into a second-order differential equation in /: 

/" - (t 2 - e)f + (a-et 2 - 2t)f = 0. (C.2) 

It is easy to show that the solution for a = reduces to the Airy integral where 
f(t) = e^/ 3 so we shall factorize / in the form: f = ug where u = e* 3//3 . The equation 
for g now simplifies tremendously: 

g" + {t 2 + e)g' + ag = 0. (C.3) 

Unfortunately there is still no closed form solution for this equation although it can 
be proved that solutions for g must be an entire function, by which we can develop an 
expansion as an infinite series: 

oo 

S = £a n r +rt . (C.4) 

n=0 

The boundary condition requires that when at — > then g — > 1. With this we can show 
that m = is the only allowed index. Then we have the following 4 term recursion 
relation: 

{n — \){na n + ea n _i) + {n — 3)a„,_ 3 + «a n _ 2 = 0. (C.5) 

The first few solutions are easily obtained: ao = 1, a 2 = f , a 3 = ^, a 4 = ^(a — e 2 ) and 
so on. As required only the first term survives in the limit a = 0. Interestingly the 4 
term relation equation (jC.5|) can actually be reduced to a three term one. Let 7„ = ^-j- 
then 

(n - l)(n 7 „ + e)7n-i + + a = 0, (C.6) 

7n-2 
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which may be more convenient for computation. Finally we have now effectively re- 
summed the series to obtain 



the last term following from the equivalence between the Airy functions and the Bessel 
functions of order 1/3. In view of the recursion properties for the derivatives of the 
Bessel functions, we have now converted the solution effectively into an infinite series of 
Bessel functions. Only the first term survives in the image free case which is now trivial 
and all the oscillatory bits are now absorbed into the Bessel functions. 




71=0 



(C.7) 



We now see that 




(C.8) 



